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ABSTRACT 

We present results from the QMAP balloon experiment, -which maps the Cosmic Microwave Back¬ 
ground (CMB) and probes its angular power spectrum on degree scales. In two separate flights, data 
were taken in six channels at two frequency bands between 26 to 46 GHz. We describe our method 
for mapmaking (removal of I//-noise and scan-synchronous offsets) and power spectrum estimation, as 
well as the results of a joint analysis of the data from both flights. This produces a 527 square degree 
map of the CMB around the North Celestial Pole, allowing a wide variety of systematic cross-checks. 
The frequency dependence of the fluctuations is consistent with CMB and inconsistent with Galactic 
foreground emission. The anisotropy is measured in three multipole bands from £ ^ 40 to £ ^ 200, and 
the angular power spectrum shows a distinct rise which is consistent with the Saskatoon results. 

Subject headings: cosmic microwave background - methods: data analysis 


1. INTRODUCTION 

The QMAP balloon experiment was designed to map 
the Cosmic Microwave Background (CMB) and measure 
the angular power spectrum on degree scales. QMAP op¬ 
erates in the Ka (^ 30 GHz) and Q-band (^ 40 GHz) with 
six detectors in two polarizations (Kal and Ka2; QI and 
Q2; Q3 and Q4), with angular resolution between 0?6 and 
0?9. Data were taken during two flights in 1996, the first 
(hereafter ELI) in June in Palestine, Texas, and the second 
(hereafter FL2) in November in Ft. Sumner, New Mex¬ 
ico. QMAP scanned a 527 square degree region near the 
North Celestial Pole in a complicated criss-cross pattern 
which allows a number of internal checks of the integrity 
of the measurements and results in an interconnectedness 
between pixels that enables efficient 1 //-noise removal. 

A detailed description of the design and performance of 
the QMAP instrument and results from the first flight are 
presented in Devlin et al. (1998, hereafter D98). QMAP 
calibrations and results from the second flight are pre¬ 
sented in Herbig et al. (1998, hereafter H98). In this Let¬ 
ter, we present the method used to analyze the QMAP 
experiment (the mapmaking process and the power spec¬ 
trum extraction), as well as the combined results from 
both flights. 

2. METHOD 

2.1. From Scan Pattern to Map 

For each of the six channels, the QMAP raw data set 
consists of M=35,318,400 observed data points yi which 
we store in an M-dimensional vector y. We subdivide the 
mapped region into N pixels whose centers form a rec¬ 
tangular grid and let Xi denote the true sky temperature 
in the direction f ^. Grouping these pixel temperatures into 


an fV-dimensional map vector x, we can write 

y = Ax-fn, (1) 

where n denotes the random instrumental noise vector of 
size M and A is a matrix of size M x N that encodes the 
QMAP scan strategy. We model n as a random variable 
with zero mean and covariance matrix given by N = (nn*). 
The scan strategy matrix A has A.y = 1 if the obser¬ 
vation points to the pixel, 0 otherwise. 

The goal is to compute a map x that estimates the true 
map X from the raw data y. We use a linear method 


X = Wy, (2) 

specified by some N x M matrix W. Substituting (^ into 
(H) shows that the error in the recovered map is 

e = X — X = [WA — I]x J- Wn, (3) 

where I is the identity matrix. Choosing W to be 
(Tegmark 1997a, hereafter T97a) 

W = [A‘M-iA]-iA‘M-i (4) 

for some MxM matrix M gives WA = I, so that x = x-f-n 


can be interpreted as an honest-to-goodness map where 
the pixel noise e = Wn is independent of x. The noise 
covariance matrix of the map is simply 

S = (ee*) = WNW*. (5) 

2.2. Pre-whitening & Pink Noise Removal 

Ideally, we would like to use the minimum-variance map¬ 
making method (Wright 1996; T97a), which corresponds 
to the choice M = N~^. In our case, the huge matrix N 
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is far from diagonal, since long-term 1// drifts (so-called 
pink noise) introduce strong correlations between the noise 
rii at different times. In other words, although direct ap¬ 
plication of equations (^) and (||) with M = would 
give what we need in principle, inverting the non-sparse 
matrix N would require a Hubble time in practice. 

Fortunately, we can obtain virtually the same answer 
by employing a series of numerical methods detailed in 
Tegmark (1997b, hereafter T97b). Since the statistical 
properties of the noise are virtually constant in time, be., 
N is almost a circulant matrix^, we replace the raw data 
y by a high-pass filtered data set 

y = Dy, (6) 

where D is another circulant matrix (a convolution filter), 
chosen so that both D and the filtered noise covariance 
matrix N = (hn*) = DND* are band-diagonal. Wright 
(1996) referred to this as “pre-whitening” and chose the 
filter so that N « I. In our case, however, N is not quite 
circulant because of the omission of ^ 600 segments of cal¬ 
ibration data (H98), but of the form N = Nc + Ng, where 
Nc is circulant and the correction Ng is extremely sparse. 
We therefore choose M = N~^, with y as the new data 
set. 

Defining A = DA, we can rewrite ( 0 ) as y=Ax -I- h, so 
equations (||) and (||) now give 

51 = [A‘D*MDA]-iA‘D‘MDy (7) 

S = [A‘D‘MDA]-i[A*D*MNM*DA][A‘D*MDA]-i 
which can be evaluated on a workstation in about 24 hours. 

2.3. Offset Removal 

In order to make maps from the data, we need to re¬ 
move instrumental offsets that are synchronous with the 
chopper position. Although we find no evidence of atmo¬ 
spheric emission or sidelobe contamination, there is ther¬ 
mal emission from the instrument at the mK-level (H98). 
We solve for this scan-synchronous component by adding 
160 “virtual pixels” X 2 (corresponding to the 160 sampling 
positions along the scan) to the map vector xi and widen¬ 
ing the scan strategy matrix A with an additional “1” in 
each row in one of 160 extra columns. This allows us to 
rewrite equation (|^ as 

y = Aixi -I- A 2 X 2 = Ax -t- n, (9) 

where Ai and A 2 are matrices of sizes MxN and M x 160, 
respectively, so x becomes a vector of dimension N + 160 
and A a matrix of size M x {N + 160). Since the scan 
strategy is so well interconnected, this method is able to 
produce an offset-free map at the cost of only a marginal 
increase in pixel noise. However, we do not wish to assume 
that the offset remains constant during the entire flight, 
as even a 10% variation in a 1 mK offset would cause an 
artificial signal comparable to the CMB. Since our offset 
is almost entirely localized in frequency to the first two 
harmonics of the scan rate (H98), we therefore combine 
the virtual pixel method with an extremely conservative 

circulant matrix is one where each row is the previous one 
diagonalized, etc.) with Fourier methods (see T97b). 


approach illustrated in Figure 1: we make D a notch filter 
which annihilates all signals at these two frequencies, as 
well as the DC (0 Hz) component. Our filtering is thus 
more of a “pre-blueing” than a pre-whitening. The price 
we pay for this conservatism is that we lose essentially all 
information about the CMB dipole, which was detected. 
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Fig. 1.— The filtering procedure is illustrated in the frequency 
domain (left) and in the time domain (right) for the Kal channel 
of flight 2. The noise power spectrum (top left) contains a 1//- 
component which causes the noise to be almost perfectly correlated 
between measurements close together in time (top right). By multi¬ 
plying the Fourier transformed signal by an appropriate filter (mid¬ 
dle left), which corresponds to applying a convolution filter to the 
signal (middle right), we obtain filtered data y that has a white 
noise power spectrum with three “notches” (lower left), correspond¬ 
ing to a block-diagonal time autocorrelation function (lower right). 

1 sample p::: 1.36 ms. Vertical units are arbitrary. 

2.4. Combining maps 

When combining two maps xi and X 2 of the same an¬ 
gular resolution into a single map x, we use the minimum- 
variance weighting 

5c= [Sri + S^i]”bsr^xi+S^ix2] . (10) 

The resulting covariance matrix for 5c is therefore 

s = [sy+ S-1]-'. (11) 

When combining maps of different resolution, the one with 
the higher resolution was first smoothed to the lower res¬ 
olution. 

shifted one notch to the right. They are easily manipulated (inverted, 
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2.5. Wiener-Filtered Map 

In addition to the x map, we also compute the Wiener 
map x„ given by (T97a) 

x™ = S[S + N]-ix, (12) 

where S=(xx*) is the CMB covariance matrix, defined as 

= E (13) 

1=1 

To avoid imprinting features on any particular scale, we 
use a flat fiducial power spectrum Ci normalized to Q = 
30/iK, which is roughly the CMB power level we find in 
the maps. We approximate the QMAP beam by a cir¬ 
cular Gaussian with FWHM values given by D98, z.e., 
Bi K, e-«"d^+i)/2 and 6 = V81n2 FWHM. The Wiener 
filtered map x^ contains the same information as x, but 
it is more useful for visual inspection since it is less noisy. 

2.6. From Map to Power Spectrum 

The signal-to-noise (S/N) method (Bond 1995; Bunn 
and Sugiyama 1995) compresses the information content 
of a CMB map into a vector z = B*x, where B is an 
N X N matrix whose i*^ column satisfies the generalized 
eigenvalue equation 

Sbi = AiSbj, (14) 

normalized so that b*Sbi = 1 and sorted by decreasing 
Xi. As described in e.g. T97b, the N numbers Zi are un¬ 
correlated, i.e., 

(z^Zj) = b-(S -f S)hj = [1 -h Xi]6ij, (15) 

and their variance (zf) has a contribution of 1 from noise 
and Xi from signal. This means that the eigenvalue Xi 
can be interpreted as a S/JV ratio for Zi, and the quanti¬ 
ties Qi (X {zf — 1) can be used as estimators of the power 
spectrum 5Ti = (.{£ since {qi) oc WiSTg for 

some window function Wi. The qi tend to probe smaller 
scales as i increases and the S/N drops. The band power 
measurements in Table 1 and Figure 3 have been com¬ 
puted by normalizing the individual qi so that their win¬ 
dow functions integrate to unity and then averaging them 
in bands with a minimum-variance weighting, to minimize 
error bars. 

3. DATA ANALYSIS 

3.1. Pipeline Tests 

We tested our data analysis pipeline by generating mock 
raw data sets that incorporate the QMAP scan strategy. 
These mock data sets were processed though the pipeline, 
recovering the original maps. When adding Monte Carlo 
white and pink noise to these mock data sets, we recov¬ 
ered maps with pixel noise consistent with the noise co- 
variance matrix S computed by the pipeline. Likewise, 
when adding scan-synchronous offsets to these mock data 
sets, we recovered the original maps as well as the offsets. 
As expected, the original maps were faithfully recovered 
even when the first harmonics of these mock offsets were 
varied slowly throughout the flight. 


We repeated our analysis for a range of pixel sizes. As 
expected, we found that as long as the pixels were smaller 
than the Shannon oversampling limit (about 2.5 times 
smaller than the FWHM), the maps x were virtually in¬ 
dependent of the pixelization. 

We made N as band-diagonal as possible using a filter 
D of band width L, then neglected the tiny elements of N 
further than L/2 from the diagonal. Tests with increasing 
L-values showed that the results converged for L ~ 150, 
so we used L = 320 in the analysis to be conservative. 

To test whether the Wiener maps were sensitive to our 
choice of power spectrum normalization, we generated 
Wiener maps for fiducial power spectra with Q=20, 30 
and 40/rK. The visual difference between the two extreme 
normalizations was minimal: the maps had the same spa¬ 
tial features in the same locations, the 20/iK map simply 
being slightly smoothed relative to the 40/xK map. 

Finally, if the beam size 6 is overestimated by 1%, the 
band power STi is overestimated by — 2] percent. 

The first term comes from the above-mentioned Gaussian 
beam correction Bi and the second from the calibration, 
which involves the beam area oc 9^. Repeating our full 
analysis with the assumed FWHM reduced by Icr (~ 3%), 
the first effect decreases the normalization of the two com¬ 
bined Ka band powers in Table 1 by 0.3% and 4%, re¬ 
spectively, whereas the second effect of course gives a 6% 
increase. 

3.2. Data Tests 

The above-mentioned scan-synchronous offset was 
around ImK (FL2) and lOmK (FLl) peak-to-peak. Al¬ 
though our notch filter technique immunized the results 
towards drifts in this offset, no such drifts were actually 
detected. 

How statistically significant is our detection of signal in 
the maps? Consider the null hypothesis that a map x con¬ 
tains merely noise, i.e., (xx*) = S. Given the alternative 
hypothesis (xx*) = S -|- S, one can show that the most 
powerful “null-buster” test for ruling out the null hypoth¬ 
esis is using the generalized y^-statistic 

2 _ x*S-iSS^ik-trS-iS 
^ ~ [2tr 

which can be interpreted as the number of “sigmas” at 
which the null noise-only hypothesis is ruled out. The re¬ 
sults of this test are given for all the individual maps in 
D98 and H98, and show that signal is detected at signifi¬ 
cance levels above 15 ct in both flights. 

Of these 11 maps (see D98 and H98), many have a sub¬ 
stantial spatial overlap. This allows a series of powerful 
consistency tests, since many potential contaminants af¬ 
fect data from different bands and polarization channels 
differently. As detailed in D98 and H98, we applied the 
same null-buster test to the difference of the various maps 
in each band where they overlap spatially, and in all cases 
found the difference maps consistent with pure noise. Our 
S/N eigenmode analysis gives the same conclusion: sig¬ 
nificant signal in the best eigenmodes of the individual 
channels, but S/N-coefficients consistent with pure noise 
in the difference maps. Thus all of the significant signal 
appears to be common to the different channels, indicating 
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that the bulk of the detected signal is due to temperature 
fluctuations on the sky. 

3.3. Foreground Contamination 

To constrain the frequency dependence of our signal, we 
repeated the null-buster test for weighted difference maps 
of the form 

X = xi - (17) 

where the map xi and the frequency vi was for the Ka- 
band and X 2 and 1/2 was for the Q-band. This placed a 2—a 
lower limit on the spectral index /3 of —1.4, which means 
that the signal cannot be dominated by foregrounds such 
as free-free emission (/3 ~ —2.15) or synchrotron radiation 
(/3 ^ —2.8). A more detailed foreground analysis, cross- 
correlating the maps with various foreground templates, 
will be presented in a separate paper (de Oliveira-Costa et 
al. 1998). 



dashed circle shows the Saskatoon coverage. 


4. RESULTS & CONCLUSIONS 

The Wiener map obtained by combining all the data 
from both flights is shown in Figure 2. The above- 
mentioned generalized y^-test shows that the signal ob¬ 
served in this figure is significant at the > 15cr-level. 

This map covers 527 square degrees, and has a substan¬ 
tial overlap with the 200 square degree Saskatoon map. 
Visual comparison of the two maps in the overlap region 
reveals striking similarities, providing further indication 


that the bulk of the detected signal is due to temperature 
fluctuations on the sky rather than systematic effects. A 
detailed statistical comparison of the QMAP and Saska¬ 
toon data sets will be presented in a future paper. 
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Fig. 3. — Angular power specwum'S*^^omoined data set shown in 
Figure 2. The results are listed in Table 1. 


The power spectrum for the total data set (FL1-I-FL2) 
is given in Figure 3 and Table 1. We see that it agrees 
well with the Saskatoon power spectrum (Netterfield et 
al. 1997), showing a rise on degree scales to power levels 
substantially above those found on large scales by COBE. 
The calibrated raw data with pointing will be made public 
after publication of this Letter. 


Table 1. — The angular power spectrum. The band pow¬ 
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respectively. Calibration errors are not included (see H98). 
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